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Inelastic back-scattering of stray light is a long-standing and fundamental problem in high- 
sensitivity interferometric measurements and a potential limitation for advanced gravitational-wave 
detectors. The emerging parasitic interferences cannot be distinguished from a scientific signal via 
conventional single readout. In this work, we propose the subtraction of inelastic back-scatter sig¬ 
nals by employing dual homodyne detection on the output light, and demonstrate it for a table-top 
Michelson interferometer. The additional readout contains solely parasitic signals and is used to 
model the scatter source. Subtraction of the scatter signal reduces the noise spectral density and 
thus improves the measurement sensitivity. Our scheme is qualitatively different from the previously 
demonstrated vetoing of scatter signals and opens a new path for improving the sensitivity of future 
gravitational-wave detectors and other back-scatter limited devices. 


Parasitic signals, arising from inelastic back-scattering 
of stray light, are a recurrent issue in high-sensitivity 
laser interferometers like gravitational-wave (GW) detec¬ 
tors. These detectors use intense laser light to measure 
differential arm length changes in a Michelson interfer¬ 
ometer topology, reaching strain sensitivities in the or¬ 
der of 10 _22 /v / Hz [T|. The scattered light problem was 
discovered right in the beginning in the first prototypes 
[2] and has been observed in all first-generation detec¬ 
tors [3]. An overview on this topic can be found in 
[4]. A sketch of a possible scatter scenario is depicted 
in Fig. [l] Shown is the advanced GW-detector topology 
that includes power recycling, (tuned) signal recycling 
and arm cavities, placed in a vacuum chamber. Imper¬ 
fect anti-reflection coatings of transmissive optics and mi¬ 
cro roughness of mirror surfaces lead to stray light. Small 
amounts are back-reflected from the surroundings, in this 
example, from the walls of the vacuum tank whose mo¬ 
tion is excited by acoustic or seismic disturbances from 
the environment. Changes in the optical path length 
cause phase modulations of the stray light, which are 
associated with changes of the optical frequency. Due to 
this inelastic character of the scattering, the recombina¬ 
tion with the interferometer mode then produces a distur¬ 
bance signal in the output light. For motions at audio fre¬ 
quencies, these scatter signals show up directly in the de¬ 
tector’s most sensitive band. Even worse are sources that 
move at lower frequencies but with large motional ampli¬ 
tudes of several wavelengths, e.g. due to micro-seismic. 
In that case, frequency up-conversion leads to broadband 
scatter shoulders that can completely cover the most in¬ 
teresting gravitational-wave signal band. Observations of 
this kind were described for example for VIRGO’s second 
science run [5]. For mitigation, the two main approaches 
are to reduce the amount of stray light as far as possible 
and to improve the isolation of potential scatter sources 
from the environment to reduce their motion. Clearly, 


technical feasibility limits both methods. In addition, fu¬ 
ture increase of the light power for shot-noise reduction 
[6] will make the detectors even more vulnerable to in¬ 
elastic back-scattering. Also, the extension towards lower 
frequencies (< 10 Hz) in third generation detectors will 
require significantly improved mitigation schemes against 
back-scattered light 0. 

In this work we propose and demonstrate a new 
approach, complementing the existing mitigation tech¬ 
niques. We collect information about the scatter signal 
independently from the scientific phase signal, using the 
orthogonal (amplitude) quadrature. Up to now, such 
additional information has been used solely for vetoing of 
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Figure 1: Scenario for the occurrence of inelastic back- 
scattering in gravitational-wave detectors. The picture 
shows the advanced detector topology with arm cavities and 
power- and signal-recycling mirrors (PRM, SRM), placed in 
vacuum. Light is scattered out of the interferometer mode, 
e.g. due to micro-roughness of mirror surfaces, and accumu¬ 
lates a time-dependent phase shift ^sc (£) due to changes of the 
optical path length. Back-scattering into the interferometer 
mode produces a disturbance signal at the output. 
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disturbance signals [HI 13, i.e. excluding corrupted mea¬ 
surement data from scientific analysis. In our approach, 
we use the data from the amplitude measurement x(t) 
to calculate a time-dependent phase-space model of the 
scattered light 'ip sc (t) and subtract its projection from 
the phase measurement data p(t). We thereby restore 
sensitivity in the affected frequency range. 

Stray light travels an extra path with respect to the 
interferometer mode and generally leads to phase as well 
as amplitude modulations of the interferometer output. 
A displacement of the back-scatter source changes the 
optical path length s(t) and results in a time-dependent 
phase shift ip sc (t) = x s W f° r scattered light. Here, 
A denotes the laser wavelength. The projections of the 
scatter signal into the amplitude and phase quadratures 
x(t) and p(t) are given by 

x sc (t) = A sin ip sc (t) and p ac (t) = A cos ip sc (t ), (1) 

where the back-scatter signal amplitude A for instance 
depends on the transfer function of the interferometer 
and on the intensity of the scattered light. The interfer¬ 
ometer’s amplitude quadrature (x(t)) readout does not 
contain any GW-signals. If the laser input light is well- 
stabilized or the interferometer has equal arm lengths 
and is operated at a dark fringe, this readout solely con¬ 
tains back-scatter signals and shot-noise. It is thus an 
unbiased monitor for the scatter signal from which the 
back-scatter signal amplitude A and the time dependent 
phase shift ip sc (t) can be extracted by fitting an analyt¬ 
ical model of the back-scatter source to the x(t) data. 
The projection into the phase quadrature (p(t)) can be 
calculated according to the right side of Eq. 0. 

In this work, we consider the case of a scatter shoul¬ 
der that is produced by back-scatter sources with large 
motional amplitudes, like already mentioned above and 
described in mg. To keep the analytical model simple, 
let us consider a single source that is moving sinusoidally 
at a constant average distance to the interferometer. This 
results in a modulated optical path length 

s(t) = s 0 + m sin(27r f m t + ip m ) (2) 

with a constant average path length so, modulation depth 
m, frequency / m and phase ip m . The signal amplitude A 
of Eq. 0 is considered to be frequency independent and 
also constant over the time measured. The Doppler shift 
of the back-scattered light is proportional to the change 
of the optical path length fds(t) = yx which results in: 

27r 

fds(t) = — m f m cos( 27r f m t + (p m ). (3) 

The maximum frequency component as observed in a 
single-sided spectrum is then given by /^ ax = | x m f m \ 
and for a modulation depth m > \/2i r, this represents 
frequency up-conversion. 


Although the model described here is quite simple, 
it already reproduces the basic structure of typical dis¬ 
turbance signals as observed in GW-detectors. It does 
generate a broadband ‘shoulder’, which even includes an 
overlaid bump structure, as observed in the LIGO data 
presented in Fig. 3.1 of reference [3]. This might indicate 
that indeed slowly varying mechanical vibrations give rise 
to a significant share of the back-scattered light distur¬ 
bances observed in GW-detectors. Our technique is not 
restricted to a simple sinusoidal model, in fact, our exper¬ 
imental demonstration actually required a model that in¬ 
cluded higher harmonics in the motion of the back-scatter 
source. 

In the experiment, we prepared a table-top Michel- 
son interferometer that was limited by inelastic back- 
scattering over a broad frequency range, as described be¬ 
low. The setup is depicted in Fig. [2] The interferometer 
had an arm length of about 7 cm and an input power of 
10 mW at a wavelength of A = 1064 nm. It was stabilized 
to a dark fringe using the Schnupp-modulation technique 
m- To artificially produce a scatter shoulder, according 
to Eq. 0 - 0 , we picked off a low power light beam in 
front of the interferometer and injected it through the 
back of end mirror M e in the east arm. We modulated 



Figure 2: Schematic of the experimental setup. In 

our table-top Michelson interferometer two signals were being 
generated: A scientific signal (pure phase modulation) in the 
north arm and a parasitic signal (broadband scatter shoul¬ 
der, appearing in phase and amplitude) in the east arm with 
the help of an electro-optic modulator (EOM) and a piezo 
actuated mirror M sc . The output signal was split at a 50/50 
beam splitter. Two balanced homodyne detectors (BHD1&2) 
measured the orthogonal quadratures (p(t),x(t)). ‘0°’refers 
to the local oscillator phase that is optimal for the scientific 
signal. DAQ: data acquisition system. 
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its phase with a piezo actuated mirror M sc , driven sinu¬ 
soidally at a frequency of 5 Hz and with a modulation 
depth of a few A , which resulted in a maximum fre¬ 
quency shift of /^ ax ~ 0.2 kHz. The center of motion 
of the mirror M sc was not stabilized but turned out to 
be sufficiently constant over the measurement time. An 
additional phase modulation at 5.2 MHz was imprinted 
by an electro-optic modulator (EOM) to shift the scatter 
signal into the MHz range. At these frequencies our mea¬ 
surements were generally limited by optical shot-noise 
and we could avoid additional disturbances, e.g. due to 
acoustics. Note, that the light was picked off outside the 
interferometer because this way it was easier to produce 
sufficiently large scatter signals in our table-top experi¬ 
ment. We also introduced a scientific signal due to a real 
differential arm length change, to show that these kind of 
signals are not affected by the subtraction of the scatter 
model. The scientific signal we generated by modulating 
the piezo mounted end mirror M n in the north arm of 
the interferometer with a sound file m- The file con¬ 
tained about 4.5 seconds of a simulated inspiral of two 
neutron stars with equal masses. The sound file was fed 
into an Agilent 33500B series waveform generator as an 
external modulation and shifted by 5.2 MHz before it was 
put on the piezo. To measure the phase and amplitude 
quadratures simultaneously, the output signal was split 
up at a 50/50 beam splitter. This resulted in 50% loss 
for the scientific signal. Generally, other splitting ratios 
can be used, compromising between loss for the scientific 
signal and signal-to-shot-noise-ratio in the scatter mon¬ 
itor. The outputs of the beam splitter were read out 
with two balanced homodyne detectors (BHD1&2) each 
using a local oscillator power of about 8 mW. The detec¬ 
tors were stabilized 90 ° out of phase and the phase space 
orientation with respect to the interferometer signal was 
adjusted by minimizing a marker peak at 5.2 MHz+1 kHz 
in the live spectrum of BHD2. The marker was also gen¬ 
erated with the piezo actuated end mirror M n . The de¬ 
tected signals were mixed down at 5.2 MHz to recover 
the audio-band signals and sent through an anti-aliasing 
filter. They were then acquired with a PCI-6259 card 
from National Instruments and processed in Lab View. 
The post-processing was done in Matlab. The measure¬ 
ment results, employing the dual readout for the scatter 
limited interferometer, are depicted in Fig. [3] 

Fig. |3](a) shows the measurements from both detec¬ 
tors in spectrograms to cope with the time dependence 
of the injected signals. The scatter signal shows up as 
‘arches’ in these plots, which are more clearly visible in 
the zoomed-in cutout of the data from BHD2 (blue). 
The arches follow from the absolute value of the time- 
dependent frequency shift of Eq. ©• The injected GW- 
signal is almost completely covered by the scatter signal, 
only its ‘tail’ towards high frequencies is clearly visible 
in the data of BHD1 (red). In the averaged spectrum of 
Fig. [3](b) the scatter signal appears as the name-giving, 


broadband shoulder in both phase (BHD1, solid red) 
and amplitude (BHD2, solid blue) quadrature. Its bump 
structure results from the projection into the quadratures 
because the different frequency components are being 
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Figure 3: Dual readout data (demodulated at 

5.2 MHz). (a) Spectrogram of the measurement data from 
BHD1V2 (red V blue) showing the time resolved scatter shoul¬ 
der (arches) concealing the simulated gravitational wave sig¬ 
nal (inspiral) that was injected from an audio file [Tl]. (b) Av¬ 
eraged spectrum showing the scatter shoulder limiting both 
quadrature measurements over a bandwidth of about 0.2 kHz. 
Apart from this, the measurements were limited by optical 
shot noise, (c) Original time data of BHD2 (averaged) with 
two fitted scatter models. A higher order model provides the 
better fit because it takes the piezo nonlinearity into account. 







































4 


generated at different total distances from the interfer¬ 
ometer. The averaged power spectral density (psd) was 
computed with Mat lab’s ’pwelch’ function, using a Han¬ 
ning window spanning half the oscillation period of the 
scatter source (At = l/(2/ m )) and an overlap of 50%. 
These settings were chosen to get a sufficient frequency 
resolution without reducing the shoulder to the harmon¬ 
ics of / m by temporal averaging. The scatter shoulder is 
clearly the dominant noise source for frequencies below 
0.2 kHz, whereby this frequency refers to a demodulated 
frequency from the MHz regime. Above 0.2 kHz the mea¬ 
surements were limited by optical shot-noise. The peak 
at 1 kHz in the spectrum of BHD1 shows the earlier men¬ 
tioned marker, used to determine the quadratures. The 
injected GW-signal can not be identified in this plot, only 
a slight rise of the noise floor at the edge of the scatter 
shoulder of BHD1 is visible. 

In the attempt to model the scatter signal measured at 
BHD2 with Eq. 0 & Q it turned out that the harmonic 
model for the source displacement in Eq. 0 was not 
sufficient. An averaged section of the recorded time data 
from BHD2 is depicted in Fig.[3](c) (solid gray), together 
with the obtained fit (dashed black). Especially around 
the turning points of the piezo at about 0.03 s and 0.13 s, 
there is a quite strong deviation between the harmonic 
model and the measurement data. To account for the 
nonlinear response of the piezo to the incoming sine wave, 
we included higher harmonics of the scatterer’s oscillation 
frequency in the model for the displacement 
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s(t) = y] m n sin(27T f m t + (4) 

n =0 

which leads to the computed fit given by the solid blue 
line in Fig. |( c). For the fitting, no prior knowledge 
about the parameters of the scatter source was assumed. 
For the projection of the scatter model to the phase 
quadrature, we allowed the signal amplitudes at the two 
detectors to differ by a constant factor, to compensate 
for example for an unbalanced splitting. Also, an addi¬ 
tional constant phase was added in the cosine of Eq. 0 
to compensate for a non-perfect quadrature orientation 
of the detectors. The two new parameters were fitted 
using the data of BHD1, while all parameters obtained 
in the fit of the amplitude quadrature were kept fixed. 
The final models for both quadratures are given by the 
dashed lines of the respective colors in Fig. [3](b). 

The results for the subtraction of the modeled inelas¬ 
tic back-scattering (in time domain) from the phase mea¬ 
surement at BHD1 are presented in Fig. [4j It is clearly 
visible in the spectrogram of Fig. [4] (a) that the injected 
GW-signal could be fully recovered. The arches are not 
visible anymore. In Fig. [4] (b) the subtracted data (solid 
red) is compared to a reference measurement (solid dark 
gray) recorded while the scattering was blocked and only 
the simulated GW-signal was being injected. For corn- 
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Figure 4: Results after subtraction of the scatter 

model. (a) Spectrogram of the measurement data from 
BHD1 after subtraction of the scatter model. The simulated 
gravitational wave signal m was clearly recovered, (b) Av¬ 
eraged spectrum showing the recovered GW-signal after sub¬ 
traction of the scatter model (solid red) in comparison with 
the original data (dashed gray) and a reference measurement 
where only the GW-signal was injected and no scattering 
(solid dark gray). 

parison, the dashed gray line in Fig. [4] (b) shows the orig¬ 
inal measurement data. In the scatter limited frequency 
range below 0.2 kHz, a sensitivity improvement of more 
than one order of magnitude was achieved, with a final 
sensitivity limited by optical shot-noise. 

To conclude, we present a proof-of-principle experi¬ 
ment, in which we show for the first time, to the best 
of our knowledge, that the noise spectral density of a 
measurement device can be reduced by using an addi¬ 
tional readout that detects the observable that is orthog¬ 
onal to the one that carries the scientific signal. Our 
result is thus of high relevance for metrology in general. 
We explicitly show that a broadband disturbance spec¬ 
trum due to parasitic inelastic back-scattering can be re¬ 
moved using the additional data. The balanced splitting 
of our dual readout increases the signal-normalized shot 
noise by a factor of 2 in power. This disadvantage can 
be mitigated by reducing the amount of power tapped 
for the scatter readout, which in turn should result in 
a reduced quality of the back-scatter subtraction. Com¬ 
promising the two effects requires further investigations. 
The increased shot noise is not an issue at all, however, if 
the overall noise is dominated - apart from back-scatter 
- by e.g. quantum back-action or thermal noise. Since 
a shot-noise limited sensitivity is usually only achieved 
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above some corner frequency, a carrier-tuned optical cav¬ 
ity could be used as a filter for optimizing our dual read¬ 
out. Sideband frequencies above the filter linewidth are 
reflected from the filter and should be analyzed by a con¬ 
ventional single readout, while the transmitted lower fre¬ 
quencies should be analyzed by the dual readout. 

In principle our scheme can directly be used in fu¬ 
ture GW detectors. The additional readout requires 
an optical local oscillator for balanced homodyne detec¬ 
tion, which has not been used in GW detectors so far, 
but was already suggested for future sensitivity improve¬ 
ments [12]. Presumably, all optics that guide the local 
oscillator and the interferometer output field need to be 
suspended to avoid differential phase fluctuations. 

Our result is not restricted to GW detectors, but can 
be seen as a fundamental approach to remove measure¬ 
ment noise that also reveals itself in the observable that 
does not commute with the actual measurement quantity. 
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